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I describe a procedure by which one can transform scattering amplitudes computed in the four 
dimensional helicity scheme into properly renormalized amplitudes in the 't Hooft-Veltman scheme. 
I describe a new renormalization program, based upon that of the dimensional reduction scheme and 
explain how to remove both finite and infrared-singular contributions of the evanescent degrees of 
freedom to the scattering amplitude. 

I. INTRODUCTION 

The Four Dimensional Helicity (FDH) scheme CD |H is widely used for computing QCD corrections at 
next-to-leading order in perturbation theory. It is particularly convenient for use with the helicity method 
and the techniques of generalized unitarity. Unfortunately, as I have recently shown Q, the FDH is not 
a unitary regularization scheme. The standard renormalization prescription [2] fails to remove all of the 
ultraviolet poles, leading to incorrect results at two loops and beyond. Thus the FDH cannot be viewed as a 
regularization scheme in which one can compute scattering amplitudes. Instead, it should be looked upon 
as a shortcut for obtaining scattering amplitudes in a unitary regularization scheme. Indeed, this is how the 
FDH has always been used at one-loop; final results have always been presented in the 't Hooft-Veltman 
(HV) scheme [4] using the prescription of Kunszt, et al. (21 to transform the FDH scheme result, but it was 
not clear whether this conversion was necessary or merely expedient, allowing one to match onto standard 
definitions of the running coupling, etc. 

It is now certain that one must convert the results of a calculation in the FDH scheme into results in a 
properly defined scheme. A first step in this direction was taken by Boughezal, et al. |[6l, who put forward 
a prescription for constructing the correct counterterms for renormalization. For inclusive calculations, 
performed using the optical theorem, like those considered in Refs. BlIfD, such a prescription is sufficient. 
Experiments, however, measure differential cross sections, and the power of the FDH scheme is that it 
facilitates the calculation of loop-level amplitudes, giving access to the differential information they contain. 
To make use of the full amplitude, one must control of both the infrared and ultraviolet structure. 
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In this paper, I will exploit the close relationship between the FDH and the dimensional reduction 
(DRED) [7] schemes to develop a prescription for transforming FDH scheme amplitudes, which may be 
easier to compute using unitarity methods, into HV scheme amplitudes that can actually be used in calcu- 
lations. The plan of the paper is: In Section [n] I will review the regularization schemes that will be used; 



in Section III I will review the infrared structure of QCD amplitudes; in Section IV I will define the FDH 
scheme in terms of the DR scheme, compute the anomalous dimensions that control the ultraviolet and 
infrared structure of DR scheme amplitudes through two loops and specify the procedure for transforming 
FDH scheme results into HV scheme amplitudes. 



II. REGULARIZATION SCHEMES 

All of the schemes that I will be working with are variations on dimensional regularization [4] , which 
specifies that loop-momenta are to treated as D m = 4 — 2 e dimensional. In dimensional regularization, 
the singularities (both ultraviolet and infrared) that appear in four-dimensional calculations are transformed 
into poles in the parameter a. The ultraviolet poles are removed through renormalization, while the infrared 
poles cancel when one performs "sufficiently inclusive" calculations. 



A. The 't Hooft-Veltman and conventional dimensional regularization schemes 

In the original dimensional regularization scheme [4], the HV scheme, observed states are treated as 
four-dimensional, while internal states (both their momenta and their spin degrees of freedom) are treated 
as D m dimensional. Internal states include states that circulate inside of loop diagrams as well as nominally 
external states that have infrared overlaps with other nominally external states. It turns out that one can treat 
internal fermions as having exactly two degrees freedom, just as they have in four dimensions, even though 
their momenta are D m dimensional, but massless internal gauge bosons must have (D m — 2) spin degrees of 
freedom, while massive internal gauge bosons have (D m — 1). 

The conventional dimensional regularization (CDR) scheme [8] is closely related to the HV scheme. 
In the CDR scheme, all states and momenta, both internal and observed, are taken to be D m dimensional. 
This often turns out to be computationally more convenient, especially in infrared sensitive theories like 
QCD, since one set of rules governs all interactions. Because the HV and CDR schemes handle ultraviolet 
singularities in the same manner, their behavior under the renormalization group, anomalous dimensions, 
running coupling, etc., are identical. 

In the HV and CDR schemes, internal momenta are taken to be D m = 4 — 2s dimensional. In general, 
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£ is a complex number and it's exact value is unimportant, but taking e to be real and positive (negative) 
is preferred by ultraviolet (infrared) power-counting arguments. It is important, however, that the D m - 
dimensional vector space in which momenta take values is larger than the standard four-dimensional space- 
time. This means that the standard four-dimensional metric tensor T]^ v spans a smaller space than the 
D m dimensional metric tensor, and the four-dimensional Dirac matrices y form a subset of the full 
7^. These considerations are of particular importance when considering chiral objects involving 75 and the 
Levi-Civita tensor, but cannot be neglected when, as in the HV scheme, one restricts observed states to be 
strictly four-dimensional. 

B. The dimensional reduction Scheme 

The DRED scheme was devised for application to supersymmetric theories. In supersymmetry, it is 
essential that the number of bosonic degrees of freedom is exactly equal to the number of fermionic degrees 
of freedom. In the DRED scheme, the continuation to D m dimensions is taken as a compactification from 
four dimensions. Thus, while space-time is taken to be four-dimensional and particles have the standard 
number of degrees of freedom, momenta are regularized dimensionally and span a D m dimensional vector 
space which is smaller than four-dimensional space-time. 

Because the Ward Identity only applies in the D m dimensional vector space in which momenta are 
defined, the extra 2 e spin degrees of freedom of gauge bosons are not protected by the Ward Identity and 
must renormalize differently than the 2 — 2 e degrees of freedom that are protected. In supersymmetric 
theories, the supersymmetry provides the missing part of the Ward Identity which demands that the 2e 
spin degrees of freedom be treated as gauge bosons. In non- supersymmetric theories, however, they must 
be considered to be distinct particles, with distinct couplings and renormalization properties. These extra 
degrees of freedom are referred to as "e-scalars" or as "evanescent" degrees of freedom. 

Since the evanescent degrees of freedom are independent of the gauge bosons, their self-couplings and 
their coupling to fermions are independent of the gauge coupling and of one another. The quartic self- 
coupling splits into multiple independent terms; if the gauge theory is SU(2), there are two independent 
quartic self-couplings, in SU (3), there are three independent quartic self-couplings, and if the gauge theory 
is SU (N);N > 4, there are four independent quartic self-couplings [9]. These new couplings run differently 
from the gauge coupling under the renormalization group and cannot consistently be identified with it. 

Notwithstanding its semantic appeal, the insistence on a proper compactification, so that D m C 4 in the 
DRED scheme, is problematic when dealing with chiral theories iPTOl . Chirality is a four-dimensional con- 
cept and one cannot consistently define chiral operators in a vector space with fewer than four dimensions. 
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One way around this is to adopt a hierarchy of vector spaces D s D D m D 4 (where D m = 4 — 2 e and D s is 
assigned the value D s = 4), as in the FDH scheme (described below). In such a scheme, chiral operators can 
be defined in the four-dimensional subspace of D m , just as they are in the HV/CDR schemes. Stockinger 
and Signer ifTTl [T2l have long advocated that this is the proper definition of the DRED scheme. Aside 
from the treatment of chiral operators, there are no important computational distinctions between D m D 4 
and D m C 4. In this paper, I will adopt the D m D 4 convention and refer to this variation of dimensional 
reduction as the DR scheme. 

C. The four dimensional helicity Scheme 

In the four-dimensional helicity scheme, one again defines a vector space of dimensionality D m D 4 
(again D m = 4 — 2s), in which loop momenta take values, and a still larger vector space D s D D m , (D s = 4), 
in which internal spin degrees of freedom take values. Note that the relative numerical values of D s , D m and 
4 are not important. What is important is that as vector spaces, D s D D m D 4. 

The FDH scheme, like the HV scheme, treats observed states as four-dimensional, except, as in inclu- 
sive calculations, where there are infrared overlaps among external states. When infrared overlaps occur, 
external states are taken to be D s dimensional. 

As in the DRED scheme, spin degrees of freedom take values in a vector space that is larger than that in 
which momenta take values. It would seem, therefore, that the same remarks regarding the Ward Identity 
and the conclusion that the D x = D s — D m dimensional components of the gauge fields and their couplings 
must be considered as distinct from the D m dimensional gauge fields and couplings would apply. 

That is not, however, how the FDH scheme has been used. All field components in the D s dimensional 
space are treated as gauge fields and no distinction is made between the couplings. The reason for doing 
this is to facilitate the use of helicity amplitudes in conjunction with unitarity methods, the idea being to 
"sew together" (four dimensional) tree-level helicity amplitudes into loop-level amplitudes. While helicity 
methods can be used in the CDR scheme [13], they are most transparently and compactly represented using 
four-dimensional external states. Thus, the FDH scheme demands that the gluons circulating through loop 
amplitudes have the same number of spin degrees of freedom as the external gluons of helicity amplitudes. 

Unfortunately, this framework fails to subtract all of the ultraviolet poles [3] and generates incorrect 
results. The evanescent couplings and degrees of freedom need to be renormalized separately from their 
gauge boson counterparts, but there is no mechanism within the FDH for doing so. The errors, however, are 
only of order &(e 1 ) in NLO calculations (which is the level at which the FDH has been used in practical 
calculations to date) and therefore do not adversely affect those results. At NNLO the errors would be of 
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order &{e q ) and at N 3 LO and beyond the errors would be singular in e. 

III. THE INFRARED STRUCTURE OF QCD AMPLITUDES 

The infrared structure of QCD amplitudes is governed by a set of anomalous dimensions which allow 
one to predict, for any amplitude, the complete infrared structure fBl [T31 . These anomalous dimensions 
are known completely, in both the massless and massive cases for one and two loop amplitudes, and their 
properties beyond the two-loop level are being actively studied |[T6ti25l . 

For a general n parton scattering process, the set of partons is labeled by f = {/i},=i... n . In the formulation 
of Refs. |fT5lfT71 . a renormalized amplitude may be factorized into three functions: the jet function ^ff, 
which describes the collinear dynamics of the external partons that participate in the collision; the soft 
function Sf, which describes soft exchanges between the external partons; and the hard-scattering function 
\Hf), which describes the short-distance scattering process, 



f t (a s (n 2 ),e) S f ( Pi ,^,a s (n 2 ),e) H f ( Pi ,^,a s (n 2 ) 



(1) 



The notation indicates that \Hf) is a vector and Sf is a matrix in color space lfl4l l26l 1271 . As with any 
factorization, there is considerable freedom to move terms about from one function to the others. It is 
convenient ||T6l [T71 to define the jet and soft functions, and Sf, so that they contain all of the infrared 
poles but only contain infrared poles, while all infrared finite terms, including those at higher-order in e, are 
absorbed into \Hf). 

A. The jet function in the HV/CDR schemes 

The jet function ^ff is found to be the product of individual jet functions Jff t for each of the external 
partons, 
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(2) 



Each individual jet function is naturally defined in terms of the anomalous dimensions of the Sudakov form 
factor lfl3Tl . 

„MS' 
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+ ... , 
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where 



y£>=2Q, y$=C i K = C i 



4 X) = \cf + \C F (8 - &) , Sf« = 2/^ - |CU &, 



^-^ 2+ 3C 3 ) + c^(^ + llc 2 -^)-c, W (^ + Ic 2 

4^ + ^(^-^-^) + Q W (H + ^) + ^ W . 



(4) 



Although % and y^-, are defined through the Sudakov form factor, they can be extracted from fixed-order 
calculations Il28l434"ll . y^i is the cusp anomalous dimension and represents a pure pole term. The Sf; anoma- 
lous dimensions contain terms at higher order in e, but I only keep terms in the expansion that contribute 
poles to In (cfli). Cp = (Af? — 1)/ (2N C ) denotes the Casimir operator of the fundamental representation of 
SU(iV c ), while Ca = N c denotes the Casimir of the adjoint representation. Nf is the number of quark flavors 
and Tf = 1/2 is the normalization of the QCD charge of the fundamental representation. £ n = YX=\ 1 A" 
represents the Riemann zeta-function of integer argument n. 

B. The soft function in the HV/CDR schemes 
The soft function is determined entirely by the soft anomalous dimension matrix IV, 
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ccdr/ e 2 „ (ll 2s p \ _i , 1 I «, MS Mi) , 1 / «.r" \ hi.' Hi: 

Sf [PU W ,CC S {H )' £ J- 1 + 2£ ^) r 5 / +^2l — ] r 5 f * r s f 








(5) 



4e 2 \ 7r 

In the color-space notation of Refs. |[14ll2"6ll2"7"1 . the soft anomalous dimension is given by |fT6l[17l 

$}-\LL*rTj*(£-), ltf-M>, (6) 



where A" = Ca (67 /18 — £2) — 10 T/Nf/9 is the same constant that relates the one- and two-loop cusp anoma- 
lous dimensions. The T,- are the color generators in the representation of parton i (multiplied by (—1) for 
incoming quarks and gluons and outgoing anti-quarks). 

IV. THE FDH SCHEME AT TWO LOOPS 



The failure of the FDH scheme as a unitary regularization scheme does not mean that it is of no value in 
computing higher-order corrections beyond the next-to-leading order. Even at NLO, the FDH scheme has 
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always been used as a means of obtaining scattering amplitudes in the HV scheme. There is no reason for 
that to change at two loops. The only difference is that one must recognize that the FDH scheme result is 
not a physical scattering amplitude, but only an intermediate step toward obtaining one. 

In formulating a prescription for converting FDH scheme amplitudes into HV scheme amplitudes, the 
first problem to address, of course, is that of renormalization. One solution to the renormalization problem, 
dubbed "dimensional reconstruction," has been proposed by Boughezal, et al. [6). The idea behind dimen- 
sional reconstruction is that if one knows the the one-loop behavior of an amplitude with arbitrary (integer) 
numbers of extra spin dimensions (momenta are always D m dimensional) then the correct two-loop ampli- 
tude can be determined from the renormalization constants at different integer spin dimensions. Note that it 
is a basic assumption of dimensional reconstruction that when one is computing a two-loop amplitude, the 
tree-level and one-loop terms that contribute via renormalization are essentially trivial, and that there is no 
appreciable cost to performing extra one-loop calculations if doing so saves effort on the two-loop piece. 
The transformations that I will develop will also subscribe to this viewpoint. 

While dimensional reconstruction is a completely valid approach to the renormalization problem of the 
FDH scheme, it does have some drawbacks. One drawback is that it appears that one must determine new 
renormalization constants for each process at each order of perturbation theory. This is quite different from 
working within a renormalizable theory, where the renormalization constants can be determined in advance 
through the study of corrections to 1PI Green functions. A more serious drawback is that dimensional 
reconstruction does not address the infrared structure of amplitudes computed in the FDH scheme. 

It is certain that the infrared structure of FDH scheme amplitudes is not equal to that of HV scheme 
amplitudes. It is also clear from optical theorem calculations 0161 that once the renormalization problem is 
fixed, one could proceed with FDH scheme calculations because the infrared overlaps will sort themselves 
out. For differential calculations, one needs to know the soft and collinear factorization properties of FDH 
scheme amplitudes in order to implement a subtraction scheme, but this has already been worked out (35]- 
[371 . The problem is that all of the FDH scheme amplitudes, real and virtual, contain errors, though the 
structure of the errors is such that, after renormalization, they cancel in the inclusive sum. Even if one 
were willing to live with such circumstances, one would still want to match onto standard definitions of 
the running coupling and would have to face the fact that parton distribution functions are only available in 
the CDR scheme. A far better choice is to transform the result to a framework like the HV scheme that is 
known to be unitary and correct and which can be easily connected to the parton distribution functions. 
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A. The connection between the FDH and DR schemes 

In order to develop a rigorous set of rules for transforming FDH amplitudes, it is necessary to define 
the FDH scheme in terms of a renormalizable scheme. One can do this by exploiting the close connection 
between the FDH and DR schemes. When formulating the QCD Lagrangians in these schemes, one starts 
with the standard Yang-Mills Lagrangian and then extends the fields into D v -dimensions. In the FDH 
scheme, one proceeds directly to the development of Feynman rules involving the D s -dimensional metric 
tensor and Dirac matrices d 0. In the DR scheme, however, one first splits the gluon field into two 
independent components, the D m -dimensional gauge field and the D Y -dimensional evanescent field Il9ll38l 
l39l . The metric tensor and Dirac matrices also decompose into orthogonal components. Those new terms in 
the Lagrangian that do not involve gauge fields are assigned new, independent couplings. The evanescent- 
quark-antiquark coupling is given the value g e {g L e =A% a e ) and the quartic evanescent boson couplings are 
given values 7],- ,1=1,2,3, where r\\ represents the quartic interaction that has the same color flow as the quartic 
gluon coupling, while 172,3 represent the non-QCD-like interactions. 

Thus, all of the DR scheme interactions are contained in those of the FDH scheme, they are simply not 
labeled by independent couplings and evanescent Lorentz structures.The only exception to this statement 
concerns the quartic evanescent boson couplings. Because the evanescent bosons are not protected by gauge 
symmetry, new quartic interactions, with new color-flows among the evanescent bosons, are generated by 
higher-order corrections which must be renormalized independently of the QCD-like quartic coupling that 
appears in the classical Lagrangian. In recognition of the fact that such terms will occur, they are usu- 
ally assigned independent couplings and added to the effective DR Lagrangian. The FDH scheme doesn't 
have such couplings, but this does not present a problem. The extra quartic terms introduced to the DR 
Lagrangian clean up the renormalization procedure, but there is no reason that the couplings assigned to 
these terms could not be chosen such that they do not contribute to a DR scattering process until radiative 
corrections to the QCD-like interactions demand that they appear. 

B. The connection between the DR and CDR schemes 

From the formulation of the Lagrangians, one can also draw a connection between the structure of the 
amplitudes in the DR and CDR schemes. In particular, the DR scheme Lagrangian contains all of the 
interactions that the CDR scheme Lagrangian does, plus a host of interactions involving the evanescent 
bosons. This means that the amplitudes in the DR scheme can be partitioned into a part that is identical to 
the CDR scheme amplitude and a part that involves the exchange of one or more evanescent bosons. One 



9 



need not consider the case of external evanescent bosons since the DR scheme renormalization program 
ensures that such terms contribute to the S-matrix at order e El|40l. The DR scheme sub-amplitude that 
involves evanescent exchanges will necessarily include a spin-sum over the evanescent degrees of freedom, 
with the result that this sub-amplitude will be weighted by a factor of D x = 2e. The only way that a term 
from the evanescent sub-amplitude can make a finite (or singular) contribution to the full amplitude is if 
it is weighted by ultraviolet or infrared poles. Thus, the full evanescent contribution to an amplitude up 
to order £° is part of the universal (ultraviolet or infrared) structure of the amplitude, and is controlled by 
anomalous dimensions. This means that the evanescent contribution to an rc-loop amplitude (that is the part 
that is different from the CDR amplitude) can be determined entirely in terms of ultraviolet counterterms, 
jet and soft functions and lower-order (0 to [n — l)-loop) hard-scattering functions. Thus, with a proper 
rearrangement of terms (the DR scheme defined below), at any order n the hard-scattering functions in the 
two schemes are related by 



Clearly, if one can draw a close connection between the FDH and DR schemes, one should be able to 
develop a prescription for the direct transformation of an amplitude computed in the FDH scheme to one 
that is computed in the HV scheme. From the above considerations, it is quite simple to state the connection. 

The four-dimensional helicity scheme is the DR scheme with two extra conditions: 

1 . External states are taken to be four dimensional. 

2. The evanescent couplings (a e and r/i ) are identified with a s . 

The first condition asserts the same distinction between the FDH and DR schemes as exists between the 
HV and CDR schemes. The restriction to four-dimensional external states does not affect the anomalous 
dimensions of the theory. The ultraviolet counterterms and the jet and soft functions are unchanged. The 
only changes are to the exact form of the finite hard-scattering matrix elements. The four-dimensional con- 
dition also forbids the appearance of external evanescent states. As mentioned before, the renormalization 
program of the DR scheme ensures that evanescent external states can only contribute to the S-matrix at 
order e or higher, so this restriction is of no consequence. 

The second condition is the one that violates unitarity and renders the FDH non-renormalizable. The 
evanescent couplings need to be renormalized differently than the QCD coupling, but there is no means of 




(7) 



C. A new definition of the FDH scheme 
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doing so once the couplings have been identified. Therefore, the FDH can only be used to compute bare 
(unrenormalized) loop amplitudes. 

In the DR scheme, on the other hand, one can determine the correct ultraviolet counterterms, and the 
infrared counterterms needed to remove the evanescent contribution, leaving the HV scheme result. By 
computing these counterterms in the DR scheme and then identifying the couplings, one obtains the coun- 
terterms needed to shift from the FDH to the HV scheme. 



D. Ultraviolet counterterms for the FDH 



When working within massless QCD, it is only necessary to renormalize the couplings. It is common 
in dimensional reduction to determine ultraviolet counterterms using modified minimal subtraction (this is 
known as the DR scheme), dropping evanescent terms, even if they contain ultraviolet poles, because the 
factor of D x renders them finite. This procedure means that the renormalized coupling in the DR scheme, 
af R differs from the standard coupling a^ s that appears in HV/CDR calculations by a finite renormal- 
ization. This finite renormalization corresponds precisely to the D x /e terms that were dropped from the 
/3-function. My goal is to remove all evanescent contributions, so I will include (D x /e) n terms in my defi- 
nitions of the j3-functions and anomalous dimensions. To distinguish it from the DR scheme, I will call this 
the DR scheme. 

Because there are so many independent couplings in the DR scheme, and because they mix under renor- 
malization, the simple j3o 1 2 ... labeling of the MS scheme is insufficient. Instead, I write, 



_» A r/DR / r/DR «DR dZ ^DR dZ fr R — 

/5 DR = 2 ^ 5_ = _( e ^_ + 0, S^DR+O! «DR 

^ dp 2 n \ n Z & dap* 7 




I fig 



ijklm 




(8) 



Similar equations yield 



a? R 




(9) 



oDR _ 2 a a e _ P ^e y oDR 

Pe M dll 2 % % . fa P ^ m 

" i,j,k,l,m 

~ A rjDR „DR - 
oDR ..2 " 'Is p 'Is V" 1 oDR 

~ i,j,k,l,m \ 

The values of the coefficients through three loops (for j8 DR and j3,P R ) are given in Appendix |a| Note that 
with the rearrangement of the evanescent contributions, the terms in j8 DR that are not proportional to D x 
are identical to the coefficients of the j3-function in the MS scheme. This indicates that the renormalized 
coupling of the DR scheme coincides with that of the MS scheme. 
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The ultraviolet counterterms for FDH amplitudes are computed as follows. First, one computes the 
lower loop amplitudes in the DR scheme and then expands the bare couplings in terms of the renormalized 
couplings using the j8 -functions of the DR scheme. Finally, the evanescent couplings are identified with the 
QCD coupling and the factors of D x are evaluated (D x = 2e). 



I-^(«*))fdh = |^(a.v,a e ,J7i)) CT 



DR 



a l; ,r] l ->a s 
D x ^2e 



(10) 



This will remove all of the ultraviolet terms, including the evanescent terms that appear to be finite because 
of the factor of D r . 



E. The infrared structure of the DR scheme 



The next step is to remove the unwanted evanescent component of the infrared structure of FDH scheme 
amplitudes. As with the ultraviolet counterterms, the terms to be removed can be identified by studying 
the structure of DR scheme amplitudes. The basic form of the infrared structure in the DR scheme is the 
same as in HV/CDR, but the anomalous dimensions receive evanescent corrections. In addition, there are 
new & anomalous dimensions that depend on the evanescent couplings. Through two-loops, the corrections 
and new anomalous dimensions depend only on the fermion-evanescent coupling, not the quartic evanes- 
cent couplings. Furthermore, because the evanescent couplings are not gauge couplings, there are no new 
counterparts to the cusp or soft anomalous dimensions, which are associated with the exchange of gauge 
bosons. 

I have determined the values of the infrared anomalous dimensions in the DR scheme by the direct 
calculation of two-loop amplitudes. I first determine the anomalous dimensions for external quarks from 
the Drell-Yan amplitude. I then obtain the anomalous dimensions for external gluons from the qq — > gy 
amplitude RT1 - I43T1 . In principle, it would be easier to extract the gluon jet function by calculating the 
amplitude for gg — > H, but the Higgs - gluon coupling is governed by a set of effective operators generated 
by integrating out the top quark. This system, involving operator mixing and higher-order corrections to 
the Wilson coefficients, has been studied to high order in the CDR scheme B4l [45] , but not in the non- 
supersymmetric DR scheme. 

The calculations of the infrared anomalous dimensions as well as the wave-function and vertex correc- 
tions used to extract the j6-functions were all calculated within the same framework. The Feynman diagrams 
were generated with QGRAF R6l and the symbolic algebra program FORM ll47l was used to implement 
the Feynman rules and perform algebraic manipulations to reduce the result to a set of Feynman integrals 
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and their coefficients. The method of Ref. 11481 was used to reduce the calculation of the vertex corrections 
to propagator integrals. The full set of Feynman integrals was reduced to master integrals using the program 
REDUZE-2 Il49l . REDUZE-2 offers significant improvements over the previous version [50] and was par- 
ticularly effective at reducing the non-planar double -box integrals that contribute to the qq — > gy amplitude. 
All of the master integrals needed for these calculations are known in the literature ll5TI457l . 
The jet function in the DR scheme takes the form, 



ln < /? 5R (o?,(M 2 ),Ofe(M 2 ),e) 
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(11) 



where the anomalous dimensions in the DR scheme are 
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(12) 



Note that the QCD coupling is a 4 MS , the same coupling used in HV/CDR calculations. Since I extract the 
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anomalous dimensions from amplitude calculations, I cannot separate the order e part of the one-loop 
anomalous dimensions, which contributes at two-loops when multiplied by a /3 -function coefficient, from 
the pure two-loop & anomalous dimensions. This merely constitutes a rearrangement of terms and does not 
affect the prediction of the infrared structure. 

The soft function changes very little in going to the DR scheme. This is because evanescent exchanges 
do not add new soft anomalous dimensions, they only add corrections to the existing terms. 



(i) 

Sf 



2 / ,TTJT\ 2 



(13) 



where K = Ca (67/18 — £2) — 10/9 Tf Nf — 2/9 D x Ca is again the same constant that relates the one- and 
two-loop cusp anomalous dimensions, this time in the DR scheme. 



F. Transforming FDH amplitudes into HV amplitudes 

I have now assembled all of the pieces needed to convert bare amplitudes computed in the FDH scheme 
into renormalized amplitudes in the HV scheme. To obtain an rc-loop amplitude in the HV scheme, one 
needs 

1 . The bare «-loop amplitude in the FDH scheme. 

2. The renormalized m-loop amplitudes (m € {0, ...,n— 1}) to order £ 2 ("- m ) in the HV scheme. 

3. The jet and soft functions to order n in the HV scheme. 

4. The renormalized m-loop amplitudes (m € {0, . . . ,n— 1}) to order £ 2 (n-m) j n ^ scheme. 

5. The jet and soft functions to order n in the DR scheme. 

Note that computing the «-loop squared amplitude to order £° already required the higher-order in e contri- 
butions to the lower-loop amplitudes in the HV scheme. The conversion procedure requires them in the DR 
scheme as well. 
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The first step is to expand Eq. ([TJ by orders of a s , 



n) 



HV 



DR 



i=0 



i=0 



(0 



HV 



DR 



(15) 



I now define the "renormalized" FDH scheme amplitude as 



FDH 



Bare 



/ FDH 



CT 



/FDH 



DR 



From this I find that 



DR 



D x ->2e 



FDH £ 



Finally, using Eq. Q, I obtain 



HV 



Bare 



/ FDH 



CT 



/ FDH 



1=1 



(0 



(0 



DR 



DR 



ae,l)i-»D(s 



(16) 



(17) 



(18) 



The infrared structure of the HV scheme amplitude can be extracted from y& n ') FDB in a similar way 
or constructed directly in terms of the lower order hard scattering matrix elements and the jet and soft 
functions. 

Let me now write out explicitly the transformation of a one-loop bare amplitude in the FDH scheme, 
involving n q quarks and anti-quarks and n g gluons, into a renormalized one-loop amplitude in the HV 
scheme. Starting with 



jew 



HV 



Bare 



/FDH 



CT 



/FDH 



+s 



(1) 



I DR 



Oe,iJl-mi 
Z>j-t2£ 



(19) 



I add in the infrared parts of the HV amplitude (note that the one-loop soft functions of the HV and DR 
scheme are identical) to obtain 

/ „ms\ /„ _l„ n — \ . 



HV 



+ 





Bare 


FDH 




Bare 


FDH 






\ « ) 


(" 



+ 



n q + ng 2 fiDR 

(i)_ 



/ HV 



(^ (1) -/ tl) )^k (0) ) HV +^) 



lA Bare («f ]« 9 +» 8 -2 

"Po 



7T 



2£ 



(20) 



HV 



h« e — 2 «„ n» 
— — c A - -C F - -C A 



JT (0) ) + me) 

HV 



The first line is just the bare one-loop amplitude with standard MS ultraviolet counterterm, while the second 
line is the finite shift, broken into ultraviolet, infrared n q and infrared n g pieces, identified by Kunszt, et 
al. (5j. Beyond one loop, the transformations are not so simple and involve the structure of the amplitudes 
in addition to the identities of the external states. 
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V. CONCLUSION 

In this paper, I have described a procedure for transforming bare loop amplitudes computed in the four 
dimensional helicity scheme into renormalized amplitudes in the 't Hooft-Veltman scheme. One of the 
simplifying features of the FDH, the treatment of the evanescent states as if they were gluons, renders the 
scheme non-renormalizable. Nevertheless, the FDH can be defined in terms of a renormalizable scheme, a 
variant of the dimensional reduction scheme. Through this connection to the DR scheme, I have shown that 
the differences between amplitudes calculated in the FDH scheme and the HV scheme (up to order e~) are 
either ultraviolet or infrared in origin and are therefore part of the universal structure of the amplitude which 
is controlled by anomalous dimensions. By computing these anomalous dimensions in the DR scheme, 
defined above, through two loops, I provide concrete formulas for the transformation of the amplitudes. 

The utility of such transformations lies in the close connection between the FDH scheme and the tech- 
niques of generalized unitarity and the helicity method. These techniques are a natural fit for the FDH 
scheme, but the results need to be transformed into a renormalizable scheme so that they can be used in 
practical calculations. With the procedures described in this paper, such transformations can be performed. 

Acknowledgments: This research was supported by the U.S. Department of Energy under Contract 
No. DE-AC02-98CH10886. 
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Appendix A: DR Scheme j3 -functions 



The non-vanishing coefficients for j8 DR through three loops are: 

nf« 11 1 1 

^ = T2 C A -- 6 N f --D x C A , 

P£ R = Y A C l - YA C AN f - \c F N f - n -D x C\ , fi* = ^D x C F N f 

™ 2857 , 1415 , 205 1 , 79 , 11 , 

^ = 3456^ " 3456^ " i# C ^ + + 3456^ + 576^ 

+ ° x \- 6912 Q + 432 Ca Nf + 2304 Ca Cf Nf ) + 13824^ C * 

«x> = 1^(1 -D,) , «0 = -^Ox(2 + D x ) , /?|| 01 = -^(l-Dj 

^20020 = 64^* + A) , ^20002 = ~~ 256 ^ x ^ ~~ ^ ' 

where I omit the last three indices if they all vanish. 
The coefficients of /^P 1 * through two loops are: 

P% = \ca ~C f - l 4 N f - \d x (C a -C F ), 0™ = 3 C F , 

P% = lcl- 5 -C A C F +Cj-^C A N f + 3 -C F N f + D x ^Cl + ^C A C F - 

+°i(yi-\^c F+ ^c F ), 

02 = -\ci + \c A C F -2& F -^C F N f +D x (-^C A C F + \c^\ , 

n nR 7 9 61 3 , 1 5 / 1 , 11 \ (A2) 

A?02100 = ~~ g ( 1 ~~ D x) i A?02010 = g ( 2 + D x) i A?02001 = 4 ( 1 — At) , 
A^O 1200 = ( 1 ~~ A) , Afo 1020 = ~~ "g" ( 2 + A) , A?o 1002 = ^ (1 ~~ A) , 

^%ioi = -^(l-A), 
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The three-loop coefficients that do not involve the quartic couplings are: 

15 51 



A 



e,04 



5 15 „ 



r 3 



5 9 „ 
4 + 4 6 



+ D X 



3 
128 



32 



C 3 -QQriV/ 



32 32 



6 ) +C 2 F N f 



21 33 
32 ~ 16 



C3 



/ 



256 



C A 



128 



r 3 

C A 



91 135 „ 

— I C?, 

64 32 b 



21 3 y 

1 £3 

128 64 b 



-CaC f /V/ 
/35 69 



37 33 
256 + 64^ 3 



C A C F 



-C F /v> 



13 . 249 r Ur 3 



47 27 
128 " 32 



■TV 



41 27 
16 + 16 ?3 



(5^+6^ 



~ C A N f 



+ Dl 



64 
29 3 
512 ~ 128 
1 3 



\64 64 



461 147 „ 

1 C3 

512 64 bJ 



r 3 



^)+C A C F N f 



f 49 
' V512 



128 



/ 43 



f 189 
' ~56 

3 



9 r 

H C3 

V 256 32 bi 



-c 3 



32 I28 6 J + ^ Cf ~ Ca ^ 



(l^t + 128 C3)+4 



1024 64 ^ 
/ 109 



oDR 

3 e,13 



r 3 

C A 



25 3 
64~4 
7 3 



85 15 

32 ~T 



11 



6C3 



+4 



V1024 64 
7 
2 



•3C3 



" I ^ 3 ) - c * c *»f (i - 1 ^ 3 ) (| - 1 ^) + e\ c ^ 

63 



+ D X 



-C 



+ C A N f 



13 3 
32 + 4 
3 



+ 7£)+QC f (1 + -£ 3 )+CaC; 



I 128 16 b 



+C 3 



13 3 r 

1 C3 

128 16 b 



-C A C F N f 



C A C F 



51 
128 



16 

33 



b)-C F N f 



5 _ 105 

64 ~ "16" 
/25 



3 



29 _ 27 

32 ~y^ 3 



Vl28 + i ?3) " CA ^ 



V 128 

3 

128 ~32 



g?3 

57 



Ca ) +C 3 



<^cl(§-^)-c ICF ( 



899 45 



167 _ 207 
256 ~T6" 
273 



1024 32 



£3) +C A C F N f ( 



Csj +c A cj 

111 r 



r 3 

C A 



4355 45 



1024 32 



C3 +C^C f 



128 32 
/21071 99 



+4 



13 45 r , 
256 - 16 ?3) + 



V 1024 16 
15 



C3 



131 

128 ' 

C F /v> 



CaC f 



I8C3 



-c 3 



85 
64 



1 15 

8 ~ 16 



27 



99 



641 

256 ~ 16 
3381 261 



■AT 



256 



16 



C F 



128 32 



1 



— CiN f+ -C A C P N f+ -ClN f 



33 



1024 



C ^ + i14^ Cf "5T2 CaCf 



„5r 3025 , 12601 . 453 „ - 129 , 475 „ „ A7 

® = "4608 C - + W56 ^ ~ 128 C ^ + 64 C ' + 2304 C ^ " C * C * N ' 



151 3 „ 

h - & 

1728 4 b 



+ 



23 
32 
643 
9216 

11 
9216 



576 



r 3 

C A 



883 
1728 
5 

13824 



C A C F 



C A C F 



256 



c A c F 



144 



C 2 A N r 



19 
864 



CaC f /V> 



(A3) 
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while the three-loop coefficients that do involve the quartic interactions are: 

— Q 1 35 

& D o R 3ioo = -^{l-lD x + 6Dl) + ^N f (1 —D x ) , 

&%oio = ^ (8- 18D X - 11D, 2 ) - 2^ Nf (2 + D x ) , 

&%ooi = ^ (2-19^ + 17^) - 4 ^ Nf (l-D x ) , 

pfLoo = -y (1 ~D X ) , P% 10 = f A {2+D x ) , pff 200 i = " (1 -D,) , 

oDR 801 , . oDR ^^^/^ n \ q DR. 50V . . 

Pe,21100 = ~~ |Q24 (1 — At) , P e ,21010 = (2 + At) , PV 2 1001 = ^^j"^ (4 ~ Ay) , 

&oo = tt^t (422-553D x + 131D?) - -D x ) , 



1024 v * y 1024' 

— 5 225 

& D o R 2020 = - 3^ (652 + 136 D x - 95D 2 X ) + —N f (2 + D x ) , 



— 1 315 

&02 = ^ (394-731D x + 337D?) - —N f (\ -D x ) , 



(A4) 



/3S R 2iio = |(2-^-D?), 

— 1 135 

#02101 = -256 (622-7730,+ 151 D?) + —N f (l-D x ) , 

P%n = -™(2-D x -Dl), 

pffmo = ^ (1 -D x ) , pf? 1020 = (2+D,) , 

Pe?fl002 = (1 ~D X ) , APfllOl = _ ^2~ ^ ~~ ' 

A%300 = - ^ (1 1 - WD X - Dl) , fc% 210 = - 1| (2 - - D?) , 
ft 1 ! 1201 = (1 1 " WD X " , ft% 120 = - g (2 - A. - D?) , 
PBim = f 2 (2-D x -Dl) , p% 02 = ^ (U-25D X + UD 2 X ) , 
&%030 = I {16 + 10D X + D 2 X ) , p% 21 = ^ (2-D x -D 2 x ) , 
«ioi2 = ~(2-D x -D 2 x ), fi% m = - 2. ( 14 - 25 D x + 1 1 D 2 X ) , 

A consistent description of j6 DR and /3 e DR through three loops only requires knowledge of the /3^ R 's 
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through one loop. These coefficients are: 

odr _ ^ _ 1 _ 9 _ 1 

Ptj i, 20000 — ~g > P?j 1,02000 — g /v /' P?7i, 10100 ~ 2 ' P»Ji,01100 _ — 2 ^ ' 

o6r _ 11 1 n odr _ 9 d odr _ 7 1 n 

P?]i, 00200 — ~g~ g x > Pnu 00110 — —■L — L'x, P?]!, 00101 — 2 ~~ 2 

oDR _ 9 5r _ 1 oi5r _ 9 odr _ 1 

Prj 2l 20000 — -jg 1 P772, 02000 — 24 ^' Prj 2 , 10010 — ^ ' Pj?2, 01010 — _ 2 /' 

dr 3 i5r 1 i5r 1 

P\) 2 ,00200 = Y5 ^ "Ac) > ^1)2,00110 = 2 ^ ~ D x) , ^ 00101 = - - ( 1 - D^) , (A5) 

P^oomo = ~~ "2 3 ^ ' A-j^oooi 1 = ~~ g ( 1 ~~ D*) > £^00002 = 7^ (1 — Ac) > 

oDR _9 oDR _ ' iv7 oDR _ o , n «DR _ ^ n 

Pj] 3 , 10001 — 2 ' Pi?3, 01001 — _ 2 /' P?]3,ooiio — z + ^-f Ptj 3 , 00101 — 2 

P > r| 3 R 00020 = - (2 + D x ) , P^OOOI 1 = ~~ ~7~ ( 2 + At) > P > rj3 R 00002 = _ 7 + To At ' 
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